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Abstract 

To understand the mass spectra of charged lepton and neutrino ^4 symmetry has been 
proposed in addition with the Standard SU{2)l x U{1)y model. We break A4 symmetry 
softly and the deviation from the tri-bimaximal mixing arises due to Zee mechanism. In the 
present work, we express two mixing angles ^13 and 623 in terms of a single model parameter 
and experimental observables, such as, mixing angle 612, mass squared differences Am|]^ and 
Am32. Using the experimental values of ^23? ^12? Am^j^ and Amgg we restrict the model 
parameter and we predict ^13. This model gives rise to ^13 ~ 11° if we allow la deviation 
of 623 and 2° deviation of 612 from their best fit values. Utilizing all those constraints, we 
explore the extent of CP violation parameter Jcp in the present model and found a value 
of Jcp ~ 2.65 x 10"'^ (for la deviation of 6*23 and 2° deviation of 612) consistent with the 
other neutrino experimental results. We have studied the mass pattern of neutrino and 
neutrinoless double beta decay (PPou) parameter |(Mj,)ee| in this model. 
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An interesting way to obtain hierarchical charged lepton mass matrix along with appropri- 
ate texture of neutrino mass matrix which can accommodate present neutrino experimental 
results, namely, solar, atmospheric and CHOOZ is through the introduction of non-abelian 
discrete A4 symmetry in a model [1]. The interplay of A4 symmetry predicts diagonal and 
hierarchical charged lepton mass matrix in addition with the neutrino masses which could be 
quasi-degenerate or hierarchical. Altarelli and Feruglio have proposed a version of A4 sym- 
metric model (AF model) [2] to obtain hierarchical charged lepton mass matrix along with 
tri-bimaximal neutrino mixing (sin 9i2 = 1/^3, sin ^23 = -I/V2, sin = 0) [3J. Although 
the model gives rise to ^13 = (It/esI = 0) which is consistent with the CHOOZ-Palo Verde 
experimental upper bound (6*13 < 12° at 3a), however, the non-zero value of jf/esl opens up 
a possibility to explore CP violation in leptonic sector which is the main goal of many future 
short and long baseline experiments. 

A prediction for non-zero Ues has been realized in a recently proposed modified AF model 
through the inclusion of three gauge singlet charged scalars due to radiative correction of the 
off diagonal elements of the neutrino mass matrix . The model successfully predicts solar 
and atmospheric neutrino mixing and mass-squared differences along with small but non- 
zero value of ^13 well below the present experimental upper bound for a reasonable choice 
of model parameters. A relationship between different mixing angles is an outcome of the 
model and the predictability of the model is also testable in future neutrino experiments. 

Aim of this paper is to generalize the assumptions made in [1]. In the present work we 
bring down the value of 612 to its best fit value 6*12 = 34.0° from the tri-bimaximal value of the 
above 612 = 35.26° which is at the la edge of experimental value. We investigate how much 
^13 can be within la variation of 623 and at the best fit value as well as 2° variation about 
the best fit value of 612. Then we shift our concentration to CP violating parameter Jcp and 
Dirac phase So and figure out their values. We also see the mass pattern of the neutrinos 
and also see the variation of jSjSoi, experimental parameter in this model. For our analysis 
we have used the best fit value of the mass squared differences, Aml^ = 8.0 x 10~^ eV^ and 
Am|2 = 2.1 X 10-3 eV^. 

For completeness, we briefly summarize here the model proposed in [1]. The lepton 
content with their representation and the Higgs content with their vevs and representation 
under SU{2)l x U{1)y x ^4 symmetry is presented in Table 1, where all fields except the 
charged scalar xt have been used in the original AF model [2]. The Yukawa interaction in 
the leptonic sector is given by 

L\f = yee'{(j)Tl)hd/A + y^,|^'{(j)Tlyhd/A + yrT'{(f)Tl)"hd/A + 
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Xaiilhjhu)/ + Xb{(j)slhjhu)/A'^ 



(1) 



where Xa, xi,, t/e, yfj_, Ur are Yukawa couplings and A is the new mass scale. After spontaneous 
breaking of the symmetry of the model, the above Lagrangian gives rise to the following mass 
terms 

£af = WT/A(?/ee''e + + y^rV) 

+ XaVl{u/A^){VeVe + 2z/^Z/^) 
+ Xbvl2Vs/'iJS^^{Ve^e + ^^.V^, 

+ Vrl'r - Ve^ii - ^iJ^T " ^^r^^e) + h-C (2) 

The charged lepton and neutrino mass matrices come out as 
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where 

me = yeVdVr/ A, = y^^vaVT/A, = yrVdVr/ A, 

a = 2xau/ A, d = 2xbVs/ A,mo = v'^/ A. (4) 

The charged lepton mass matrix is diagonal so the leptonic mixing solely occurs from the 
neutrino sector and diagonalising mj, by the way Wm^U* we get the three mass eigenvalues 
as 

mi = a + d, 7712 = CL, = d — a (5) 

with the exact tri-bimaximal mixing pattern. Thus, in the AF model, tri-bimaximal mixing 
occurs naturally. In |4j, a modified version of the above model has been investigated to 
generate non-zero ^13 which at the leading order gives tri-bimaximal mixing form. In order to 
do that, an ^4 triplet SU (2)^ singlet charged scalar xti= Xij xt 1 xt) has been incorporated 
and the leptonic part of the Lagrangian becomes 

L = L\f + L\iaf (6) 

where L^maf has two parts as 

L^MAF = -Ci + £2 (7) 
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Ci is symmetry preserving part and is given by 



£1 = / (L L X*) C (3 X 3 X 3) 

however, £2 contains Zee-type term which is exphcit soft A4 symmetry breaking and is given 
by 



A = Ci2hliT2hd{xt + Xt + Xt)- 



(9) 



The charged lepton mass matrix is still diagonal in the present model and the neutrino mass 
matrix comes out as 
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/a + 2d/3 -d/3 -d/3-e 
~d/3 2d/3 a-d/3 + e 

y-d/3-e a-d/3 + e 2d/3 

where the a and d parameters defined earlier and are obtained due to higher dimensional 
operators in the same way as AF model. The e term is the additional contribution arises at 
the one- loop level due to well known Zee mechanism and is shown in Fig.l. The parameter 
e is given by 
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167r2(Mf - M|) 

Although the Lagrangian given in Eq|8] can generate corrections to all off-diagonal entries of 
the mass matrix given in Eq. (fTOl) . however, dominant terms proportional to are retained. 
In this model neutrinos are Majorana-type in nature. In general, the parameters a, d, e are 
all complex, however, it is possible to rotate out one of the phase. For our analysis, we 
consider only the parameter d is complex and parameters a,e are real. The neutrino mass 
matrix in this case comes out as 



/ a + 2de'^ /3 -de''^ /3 -de'''' /3 - t 

-de''^/3 2de'^l3 a-de'^l3 + t 

y -de''^ /3 - e a- de"^ /3 + e 2de'^ /3 

A relationship between the parameters a and d has been considered as 
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Figure 1: One loop radiative f^.^i - mass due to charged Higgs exchange. 
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Table 1: List of fernaion and scalar fields used in this model. 
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where k, d and a are the real parameters. The analysis of Ref . [1] has been done for a specific 
value of K and in the present work, we have generalized the whole analysis. We constrain 
the value of n first from the existing bounds of two mixing angles 612 and 623- Then we have 
utilized the result to calculate 613 and the CP violation parameter Jqp |5j and explore the 
extent at which Jcp is allowed in the present model. 

With only assumption e is small, we diagonalize mass matrix Eq. (fT3l) as 



U^M^U* = diag (^de'^ + a + e, a, de'^ - 
upto first order in e where diagonalizing matrix 
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(16) 



It is to be noted that for vanishing value of e the matrix U leads to tri-bimaximal form. 
Three approximate mass eigenvalues in Eq. f|T5l) take the following forms 
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where e' = e/a and we use the relation given in Eq. f|T^ . The solar and atmospheric neutrino 
mass squared differences are coming out as 
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From mixing matrix U in Eq. ([TB|1 we obtain the mixing angles: 
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Figure 2: Plot of cos^0 with respect to /t and 9i2. We keep Am|2 and Aml^ to their best 
fit values. 
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We define the ratio R as 



Am| _ K cos^ (^(k + 2) + 2e'(l + k cos^ 0) 



^'^atm cos^ 0(2 — k) — 2e'(l — K cos^ < 



Using the first relation of Eq. ( fT9l) we have 

, V3/t cos^0(/t + 2) ^ / . „ ^ , 
e = — — X (sin6'i2 - l/v3) . 



2 + K COS^ y 

Using Eq. ^ in Eq. ([20D we get 

-2 r 2(1 - i?) + k(1 + i?) + V3(k + 2)(sinei2 - 1/V3)(1 + R) 



cos 



(19) 
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K [2(1 -R) + k{1 + R) + 2VS{k + 2)(sinei2 - 1/V3)(1 - R)_ ' 
Again using expression of e' from Eq. fl2T]) to the expression of Am32 in Eq. (fT8|) we get the 
dependence of the on n and cos^ 0: 

2 Am32 (2 + K COS^ 0) 



K cos2 0[(fi; - 2) (2 + K cos2 0) + 2v^(k + 2) (1 - k cos^ 0) (sin ^12 - 1 / v^)] 
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Figure 3: Plot of with respect to k and 612. We keep Am|2 and ^m^^ to their best fit 
values. 

Thus from Eq. (!22l) and Eq. (!23|) . we see that 0? only depends on single model parameter n. 
From the above expression of cos^ we can put bound on k for the given values of R and 
9i2- For the best fit value of Am|2 and Am^i we have the approximate bound: 



in the course of variation 32° < 612 < 36° for cos^ < 1 and it is shown in Fig. [21 It also 
ensures cos^ > 0. Again compatibility of the above bound of n with the restriction > 
demands that Amgg > 0. This leads to normal ordering of neutrino masses. 

Using the expression of e' from Eq. (12T|) into the expression of sin 613 and tan^ 623 in Eq. 
(fT9|) we get the two mixing angles in terms of k, cos^ 0: 
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Vq{k -2){2 + K COS^ I 



^ sin^i2-l/v^ (26) 
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Here also tan^ 6*23 and sin 6*13 only depend on k as cos^ <^ is only function of k. In addition to 
this dependence on parameter k, they also depend on the well measured quantities, mixing 
angle 612-, and the ratio of solar and atmospheric mass-squared differences R. We keep 
the solar and atmospheric mass-squared differences on their best fit values. We study the 
variation of 6*23 with mixing angle 9i2 and the parameter k in Fig|31 We have varied k 
from its analytical upper bound —2.0 to —2.5 and also have varied 612 from 32° to 36° (2° 
deviation about best fit value 34° of 612). For a fixed value of n variation of 6'23 with 9i2 is 
small and is smaller in higher value of k, e.g. for the variation of 9i2 from 32° to 36° for 
K = —2.0, ^23 remains almost at 45.8° and for k = —2.5, ^23 changes from 48.4° to 48.6°. If 
we allow ^23 upto its 1° deviated value 46°, range of k, shrinks to —2.04 < k < —2.0 which 
is insensitive to the variation of 612- To keep 6'23 within la deviated value 48° for the whole 
range of variation of 612 from 32° to 36°, we have to keep k > —2.39. From Figj3]we can 
also study the variation of 6'23 with k for a fixed value of 612, e.g. for the best fit value of 
612 = 34°, 6'23 changes from 45.8° to 48.5° for the variation of k from —2.0 to —2.5. So, from 
Figl3]we have gathered the information where we should keep the value of k in the light of 
experimental values of 6'23 and 612- With this information we study other observables. 

FigHl we have plotted 613 with 6*12 and k in their allowed region using the best fit value 
of solar and atmospheric mass squared differences. We have promised to generate nonzero 
6'i3 with changing 612 from its tri-bimaximal value 35.26°. Thus, we obtain the prediction 
on 9n as 2.91° < ^13 < 3.7° for -2.04 < k < -2.0 (alternatively for 45.8° < ^23 < 46.0°) 
for the whole range of variation of 612 32° to 36°. If we allow ^23 upto its la deviated value 
45.8° < 6*23 < 48.0° (equivalently -2.39 < k < -2.0), we have the prediction of 6^13 as 
2.91° < 613 < 10.7° for the whole range of variation of 612 from 32° to 36°. So, the upper 
bound is near the largest possible allowed value from the experiment (CHOOZ ^13 < 12° 
at 3a). From the plot of FigJH this is to be noted that the upper bound of ^^13 is mildly 
sensitive to 612, it varies from 9.93° to 10.7° at k = —2.39. But the lower bound of 613 is 
2.91° which is insensitive to the variation of 612 at k = —2.0. 

Keeping all those constraints in view next we explore the parameter space of CP violation 
parameter Jcp. The parameter Jcp defined as [5] 

J ^ ■ on ■ oa ■ OQ Q ■ X Im[hi2h23h3i] , . 

Jcp = o 26'i2 sm 26^23 sm 26'i3 cos 9i3 sm do = -r — — n"-: — n- (27) 
8 Ami-, AmiT Amin 
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where h = M^M^, Sd is Dirac phase. This Jcp is associated with CP violation in neutrino 
oscillation and is directly related to Dirac phase of mixing matrix. Using Eq. (fT3|) . Eq. (fT4l) 
in Eq. (|271) CP violating parameter Jcp takes the following form 



a 
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2e' 



upto first order term in e'. Using the expression of e' from Eq. (1211) and from Eq. 
into the Eq. fl28p we have 



2^3 + 2)2(2 + K COS^ 0) 2 (sin ^12 — 1 / \/3) \/ COS^ — COS^ 

~ 9 i?[(/t - 2) (2 + K cos2 0) + 2 v^(/t + 2) (1 - K cos^ 0) (sin - 1/^3)]^ ' 

Jcp is also only function of k, because cos^ in Eq. fl22l) is function of k only. It is clear 
from the above expression if k = —2 the value of Jcp becomes zero. This k + 2 factor appear 
in Jcp in Eq. fl28l) is purely from symmetric structure of the neutrino mass matrix. In 
our analysis we have seen that k and cos^ are not independent, k = —2 approximately 
corresponds to cos^ = 1 where Jcp also vanishes. Using the expression of cos^ from 
Eq. ( 122!) into above Eq. (!29|) we can make Jcp only k dependent. It also depends on the 
observables quantities 6*12, Amgg and Aml^. We have plotted in Figl5] Jcp with respect k 
and 612 in their allowed region keeping A'm'^2 ^"^^ ^^21 their best fit values. We have the 
values of Jcp in the range 3.27 x 10"^ < Jcp < 7.6 x 10"^ for -2.04 < k< -2.0 (alternatively 
for 45.8° < 6*23 < 46.0°) for the whole range of variation of 612 from 32° to 36°. If we allow 
6123 upto its la deviated value 45.8° < 6*23 < 48.0° (equivalently -2.39 < k< -2.0), we have 
the prediction of Jcp as 3.27 x 10^^ < Jcp < 2.65 x 10"'^ for the whole range of variation 
of 612 from 32° to 36°. Upper bound of Jcp is varying from 2.12 x 10~^ to 2.65 x 10^^ for 
the variation of 612 from 32° to 36°. For the best fit value of 612 (~ 34°), upper bound is 
2.42 X 10"'^. So, larger value of Jcp ~ 2.65 x 10^'^ is possible and it can be probed through 
upcoming base-line experiments. From the Eq. fl27j) we can find the expression for sin 5/). 
Using expressions for 6^13 from Eq. (1261) , 6*23 from Eq. (I25l) , Jcp from Eq. (l29l) and cos^ from 
Eq. (I22I) into Eq. (1271) we can have k dependent function for sin^^). We have plotted 6d with 
respect to k, in Figl6]for the best fit values of 612, Am|2 and Amji. The figure reflects the 
prediction of 6d as 6r) = 3.6° for 623 = 48° at k = —2.39. 

Now we are going to see the behavior of mass eigenvalues and their sum with respect n 
and hence the other observables. Using the the expressions for from Eq. (p3l) . e' from Eq. 
( I2TI) and cos^ from Eq. (1221) into Eq. (|T71) we get the k, dependent functions for mass eigen 
values. We have plotted mi, 1712, and their sum in FigJ7]with respect to k for the best 
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Figure 5: Plot of Jcp with respect to k and 6*12. We keep Am^g and Aml^ to their best fit 
values. 
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Figure 6: Plot of Dirac phase 5d with respect to k for the best fit values of ^12, Am|2 and 
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fit values of 6*12, Amgg and Am2i. Tlie observations of the plots in Fig 17] suggest that mass 
pattern is normal-hierarchical. We also have seen that 0.07 eV < nii + m2 + < 0.076 eV 
in the course of variation of k —2.39 < k < —2 (48° > ^23 > 45.8°). It also satisfy the 
cosmological bound mi + m2 + < 0.7 eV [6]. 

Again parameter responsible for the /3/3oj/ experiment is also studied in the present model 
and the relevant quantity: 



a + 



3 



^ 4k(k + 3) cos^ I 
^ 9 



1/2 

(30) 



Using the expressions for a and cos^ (j) from Eq. ( l23ll and Eq. ( |22l) respectively, we get |(M;^)ee| 
in terms of model parameter k and other physical observables. Keeping experimental values 
of the 612, Am32 and Amli to their best fit value, we have plotted |(Mi,)ee| with respect to n. 
For the physical region of k, FigJT] shows that |(Mj,)ee| is well below the experimental bound 
0.89 eV. 

A point to be noted as our analysis of mixing angles are matching with [4] in the real 
liniit((/) = 0°). However, the expressions for the mixing angles for complex case in are not 
in exact correspondence with our result in the present work for k = —2 limit. This is because 
in the complex analysis in [1], we have assumed that eigenvectors of M,^ construct U. This 
U has ability to diagonalize Mi^ to its true eigenvalues but this U may not be unitary. It is 
better to diagonalize h (M^Ml) which is hermitian and its diagonalizing matrix is unitary. 
But it is always not easy to do the same. In this paper we have solved the 18 equation from 
Eq. (ITSl) and find out U in Eq. (fT6!) which also is unitary keeping terms upto the order e: 
Wu = UW = 1 + 0(e2). 

In summary, we explore the parameter space of a softly broken symmetric model for 
different mixing angles and a model parameter k. We expressed the two mixing angles 612 
and 623 in terms of a single parameter n, and constrained the parameter space for the best 
fit values of Amgg and Am'^i. With the allowed parameter value we predict 61^, 61^ ~ 11° 
(for la deviation of 623 and 2° deviation of 612 about their best fit value). Utilising the 
above result, we expressed the CP violation parameter Jqp in terms of n and sin6'i2 and 
explore the extent of Jcp allowed in the present model. A comparatively larger value of Jqp 
is allowed by the present model (Jcp = 2.65 x 10~^ for la variation of the angle 623) and 
consistent with other neutrino experimental results. 
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Figure 7: Plot of neutrino masses mi, m.2, and their sum, and also \ {Mi,)^ 
to K for the best fit values of 9i2, Am|2 and Aml^. 
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